1 About

This report is based on a model presented at the 2021 International Conference on Evolving Cities, University of Southampton, 22 – 24 September 2021.

If you are interested in how the model works start from https://dataknut.github.io/localCarbonTaxModels/

The model is similar in some respects to Evans et al (2023) Getting to net zero: Islington’s social housing stock (https://journal-buildingscities.org/articles/10.5334/bc.349) which estimated an overall cost of upgrading Islington’s social housing stock of 4500 buildings containing some 33,300 dwellings to be £1,600m.

1.1 Published versions

A paper based on this model as applied to the City of Southampton has been published in Buildings & Cities:

  • Anderson, B. (2023). A residential emissions-based carbon levy: city and neighbourhood consequences. Buildings and Cities, 4(1), pp. 1–20. DOI: https://doi.org/10.5334/bc.279

1.2 Citing this report

If you wish to re-use material from this report please cite it as:

Ben Anderson (2023) Simulating a local emissions levy to fund local energy efficiency retrofit: East Hertfordshire as a case study. University of Southampton, United Kingdom

License: CC-BY

Share, adapt, give attribution.

2 Highlights

This report estimates the value of an emissions levy for East Hertfordshire using LSOA level data on emissions derived from the CREDS place-based emissions calculator. These emissions are all consumption, gas and electricity. It does this under two scenarios - a simple carbon value multiplier and a rising block tariff.

It then compares these with estimates of the cost of retrofitting EPC band dwellings D-E and F-G in each LSOA and for the whole area under study.

Key results:

Section 10 gives the overall summary of the annual £ emissions levy value and estimated retrofit costs

  • Figure 5.2 maps the LSOA level per dwelling ‘all emissions’ in Tonnes CO2e - which areas emit the most?
  • Table 6.1 shows total levy generated under Scenario 1
  • Table 6.8 shows total levy generated under Scenario 2
  • Figure 6.11 compares the scenarios in terms of % of levy generated by areas in each IMD decile while Figure 6.12 compares the levy generated under each scenario at LSOA level. In both cases, Scenario 2 should be lower in more deprived areas and higher in less deprived areas. Figure 6.11 shows whether or not this is the case
  • Table 7.1 shows total retrofit costs and Figure 7.1 shows the LSOA level retrofit costs per dwelling by IMD decile for comparison with Figure 5.1
  • Figure 9.1 shows the years to pay back under Scenario 1 for an all emissions levy while Figure 9.7 does the same for Scenario 2
  • Figure ?? shows what would happen after year 1 if the levy were shared equally across LSOAs (all emissions, Scenario 1) and Figure ?? shows the same for Scenario 2.
  • Figure 9.2 maps the years to pay under Scenario 1
  • Figure 9.12 shows payback years under each Scenario assuming a constant all emissions levy

3 Introduction

This report estimates a model of an emissions levy for East Hertfordshire using LSOA level data on emissions derived from the CREDS place-based emissions calculator.

The model applies carbon ‘values’ to a number of emissions categories to estimate the levy revenue that would be generated for each LSOA in year 1 of such a levy. It then sums these values to given an overall levy revenue estimate for the area in the case study.

The report then use estimates of the cost of retrofitting EPC band dwellings D-E and F-G together with estimates of the number of such dwellings in each of the LSOAs to calculate the likely cost of such upgrades in each LSOA and for the whole area in the case study.

Finally the report compares the distributions of the two to understand whether sufficient revenue would be generated within each LSOA to enable the per-LSOA or whole case study area costs of the energy efficiency upgrades to be met. In doing so the report also analyses the extent to which redistribution of revenue from high emissions areas (households) would be required.

It should be noted that this area level analysis uses mean emissions per household. It will therefore not capture within-LSOA heterogeneity in emissions and so will almost certainly underestimate the range of the household level emissions levy values that might be expected.

4 Data

The model uses a number of Lower Layer Super Output Area (LSOA) level datasets to analyse the patterns of emissions.

All analysis is at LSOA level. Cautions on inference from area level data apply.

## Filter: East Hertfordshire
## 
##     0     1 
## 32772    81
## Number of selected LSOAs - merged data:
## [1] 81
##    
##     East Hertfordshire
##   1                 81
## Test linkage of Census 2021 data - suspect some LSOAs have changed
##    c2021_nhhs     c2021_pc_elec_ch_2021    IMDScore       pc_fuelPoor   
##  Min.   : 478.0   Min.   : 1.541        Min.   : 1.501   Min.   : 3.00  
##  1st Qu.: 593.0   1st Qu.: 3.325        1st Qu.: 3.648   1st Qu.: 8.00  
##  Median : 672.5   Median : 7.795        Median : 7.182   Median :10.00  
##  Mean   : 704.0   Mean   : 8.775        Mean   : 8.158   Mean   :10.37  
##  3rd Qu.: 802.0   3rd Qu.:10.661        3rd Qu.:10.724   3rd Qu.:13.00  
##  Max.   :1083.0   Max.   :45.506        Max.   :27.148   Max.   :19.00  
##  NA's   :5        NA's   :5
## # NAs are LSOAs that didn't match - only matters to analysis using Census 2021 data on electricity central heating

4.1 Case study context

This report takes East Hertfordshire as a case study.

  • Number of households (Census 2021): 53,505
  • Number of electricity meters (2018): 61,397
  • Number of gas meters (2018): 50,756

Figure 4.1 shows the distribution of deprivation across the case study area using the IMD 2019 data.

If there are gaps in these maps then there may be LSOA level data linking errors, possibly due to name or boundary changes. We love name and boundary changes. If so the model results may be partial. #YMMV

Figure 4.1: IMD (LSOAs)

Figure 4.2 shows the distribution of fuel poverty across the case study area using the BEIS/DESNZ 2019 data.

Figure 4.2: Fuel poverty (LSOAs)

Observed EPCs - need to check bias with respect to e.g. Census data. EPCs more likely to be from new houses, rentals, re-sales (required) etc.

## Saving 7 x 5 in image
## Number of EPC coded dwellings in case study area, (East Hertfordshire) = 34132

Imputed EPCs (see logic in Section 7.1) - this corrects for numbers but not the bias in the EPC sample

## Number of imputed EPC coded dwellings in case study area, (East Hertfordshire) = 61389
##    case_study mean_MSOA_totInc
## 1:          0         43881.51
## 2:          1         51861.73

## Saving 7 x 5 in image

## Saving 7 x 5 in image

5 LSOA level emissions estimates

See https://www.creds.ac.uk/why-we-built-a-place-based-carbon-calculator/

“The highest carbon areas have an average per person footprint more than eight times larger than the lowest carbon areas.”

“We are not effectively targeting decarbonisation policies in high carbon areas. For example, the recently collapsed Green Homes Grants scheme provided a grant to cover 66% of the cost (up to £5,000) of retrofitting homes. For people claiming certain benefits, the cap was raised to 100% and £10,000. But the calculator shows that the big polluters are the large homes in very wealthy areas. In these neighbourhoods, the issue is not affordability but motivation. For high income households, energy costs are a small proportion of their expenditure and so the cost savings for retrofitting their home are inconsequential. As there are no policy “sticks” to incentivise action in the collective interest it is unsurprising that high carbon neighbourhoods have not prioritised decarbonisation.”

Source: https://www.carbon.place/

Notes:

  • Emissions are presented as per capita…
  • Appears to be based on residential/citizen emissions only - does not appear to include commercial/manufacturing/land use etc
Table 5.1: Summary of CREDS per capita data for LSOAs: East Hertfordshire

region

nLSOAs

mean_KgCo2ePerCap

sd_KgCo2ePerCap

East

81

10,071.5

3,104.4

## Do we have any missing energy data?
Table 5.1: caption

LAD11NM

LSOA11CD

LSOA01NM

nElecMeters

nGasMeters

CREDS_total_kgco2e

CREDS_elec_kgco2e2018

CREDS_gas_kgco2e2018

CREDS_otherEnergy_kgco2e2011

CREDS_allHomeEnergy_kgco2e2011

Now we need to convert the per capita to totals and then use the number of electricity meters as a proxy for the number of dwellings

Ideally we’d have Census 2021 data but we don’t have it yet. So instead we’ll use the number of electricity meters for 2018 which aligns with the CREDS data (might be an over-estimate if a dwelling has 2…)

First check the n electricity meters logic…

##                   LSOA11NM                     WD18NM nGasMeters nElecMeters epc_total
## 1: East Hertfordshire 001B                Buntingford       1525        1439      1050
## 2: East Hertfordshire 007D Bishop's Stortford Central       1077        1334      1040
## 3: East Hertfordshire 006A                 Puckeridge        991        1252       667
## 4: East Hertfordshire 017D            Hertford Castle        970        1676      1360
## 5: East Hertfordshire 016B         Hertford Kingsmead        964        1004       605
## 6: East Hertfordshire 004C   Bishop's Stortford Meads        940        1027       586
##                   LSOA11NM                     WD18NM nGasMeters nElecMeters epc_total
## 1: East Hertfordshire 017D            Hertford Castle        970        1676      1360
## 2: East Hertfordshire 001B                Buntingford       1525        1439      1050
## 3: East Hertfordshire 007D Bishop's Stortford Central       1077        1334      1040
## 4: East Hertfordshire 006A                 Puckeridge        991        1252       667
## 5: East Hertfordshire 004B   Bishop's Stortford Meads        773        1234       984
## 6: East Hertfordshire 010E            Watton-at-Stone        918        1118       563
Table 5.2: LSOAs with highest number of gas meters (after cleaning)

LSOA11NM

WD18NM

nGasMeters

nElecMeters

epc_total

East Hertfordshire 001B

Buntingford

1,525

1,439

1,050

East Hertfordshire 007D

Bishop's Stortford Central

1,077

1,334

1,040

East Hertfordshire 006A

Puckeridge

991

1,252

667

East Hertfordshire 017D

Hertford Castle

970

1,676

1,360

East Hertfordshire 016B

Hertford Kingsmead

964

1,004

605

East Hertfordshire 004C

Bishop's Stortford Meads

940

1,027

586

Check that the number of electricity meters reasonably correlates with the number of EPCs from the CREDS data. We would not expect the number of gas meters to correlate due to non-gas dwellings etc.

There may also be difficulties where there are multiple meters per property - e.g. one ‘standard’ and one ‘economy 7’. Really should switch to using address counts from postcode file.

Check that the assumption seems sensible…

Check for outliers - what might this indicate?

5.1 Estimate per dwelling emissions

We want to present the analysis in ‘per dwelling’ or ‘per household’ terms so we need to convert the total kg CO2e values to per dwelling values by dividing by the number of electricity meters.

## # Summary of per dwelling values
Table 5.3: Data summary
Name …[]
Number of rows 81
Number of columns 9
Key NULL
_______________________
Column type frequency:
numeric 9
________________________
Group variables None

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
CREDS_total_kgco2e_pdw 0 1 24463.55 9043.14 9915.70 17653.55 24741.09 28453.85 48240.00 ▆▇▇▂▂
CREDS_gas_kgco2e2018_pdw 0 1 2599.28 699.17 756.66 2214.85 2595.92 3000.87 4359.38 ▁▃▇▃▂
CREDS_elec_kgco2e2018_pdw 0 1 1230.44 308.44 789.58 1034.22 1123.83 1367.99 2309.55 ▇▇▃▁▁
CREDS_measuredHomeEnergy_kgco2e2018_pdw 0 1 3829.72 708.74 2386.07 3426.46 3758.43 4175.43 5886.56 ▃▇▆▃▁
CREDS_otherEnergy_kgco2e2011_pdw 0 1 234.05 409.96 0.00 33.05 57.65 152.86 1795.75 ▇▁▁▁▁
CREDS_allHomeEnergy_kgco2e2018_pdw 0 1 4063.77 789.34 2422.17 3519.25 3972.97 4714.80 5959.12 ▂▇▇▆▂
CREDS_car_kgco2e2018_pdw 0 1 3039.31 766.86 1044.19 2516.41 2981.93 3589.46 4843.42 ▂▆▇▅▂
CREDS_van_kgco2e2018_pdw 0 1 423.13 260.69 147.51 250.51 343.97 525.59 1270.27 ▇▂▂▁▁
CREDS_personalTransport_kgco2e2018_pdw 0 1 3462.44 893.10 1651.91 2818.96 3377.18 4108.11 5552.27 ▂▇▆▅▂

Examine patterns of per dwelling emissions for sense.

5.1.1 All emissions

Figure 5.1 shows the LSOA level per dwelling ‘all emissions’ in Tonnes CO2e as estimated by the CREDS tool against the Index of Multiple Deprivation (IMD) score and uses the size of the points to represent the % of dwellings with electric heating. Colour is used to represent the IMD decile where decile 1 is the 10% least deprived.

## Per dwelling T CO2e - all emissions
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Scatter of LSOA level all consumption emissions per dwelling against IMD score

Figure 5.1: Scatter of LSOA level all consumption emissions per dwelling against IMD score

## Correlation test (pwcorr)
## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$IMDScore and selectedLsoasDT$CREDS_total_kgco2e_pdw
## t = -4.5586, df = 79, p-value = 1.853e-05
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.6135674 -0.2643489
## sample estimates:
##        cor 
## -0.4563572
## Total emissions per dwelling (LSOA level) summary
##    LSOA11CD            WD18NM                         IMD_Decile_label All_Tco2e_per_dw
##  Length:81          Length:81          10 (10% least deprived):36      Min.   : 9.916  
##  Class :character   Class :character   9                      :18      1st Qu.:17.654  
##  Mode  :character   Mode  :character   6                      : 8      Median :24.741  
##                                        8                      : 8      Mean   :24.464  
##                                        7                      : 7      3rd Qu.:28.454  
##                                        4                      : 3      Max.   :48.240  
##                                        (Other)                : 1
Table 5.4: Highest emitting LSOAs (per dwelling)

LSOA11CD

WD18NM

IMD_Decile_label

All_Tco2e_per_dw

E01023452

Bishop's Stortford Silverleys

10 (10% least deprived)

48.2

E01023490

Hertford Sele

9

47.2

E01023468

Great Amwell

10 (10% least deprived)

43.4

E01023461

Bishop's Stortford South

9

43.4

E01023471

Hertford Bengeo

10 (10% least deprived)

41.9

E01023481

Hertford Heath

8

40.2

Table 5.4: Lowest emitting LSOAs (per dwelling)

LSOA11CD

WD18NM

IMD_Decile_label

All_Tco2e_per_dw

E01023487

Hertford Sele

4

9.9

E01023450

Bishop's Stortford Meads

7

10.8

E01023499

Sawbridgeworth

6

10.8

E01023488

Hertford Sele

4

11.0

E01023478

Hertford Castle

9

11.5

E01023446

Bishop's Stortford Central

8

11.5

We should not be surprised that emissions are negatively correlated with deprivation. If you are, you should try reading an excellent paper by Milena Buchs and Sylke Schnepf, both formerly of the University of Southampton:

“whilst all types of emissions rise with income, low income, workless and elderly households are more likely to have high emissions from home energy than from other domains, suggesting that they may be less affected by carbon taxes on transport or total emissions. This demonstrates that fairness implications related to mitigation policies need to be examined for separate emission domains.”

With that in mind, read on.

Figure 5.2: Annual mean T CO2e per dwelling (LSOAs)

5.1.2 Home energy use

Figure 5.3 uses the same plotting method to show emissions per dwelling due to gas use.

## Per dwelling T CO2e - gas emissions
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   756.7  2214.8  2595.9  2599.3  3000.9  4359.4
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Scatter of LSOA level gas per dwelling emissions against IMD score

Figure 5.3: Scatter of LSOA level gas per dwelling emissions against IMD score

## Correlation test (pwcorr)
## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$IMDScore and selectedLsoasDT$CREDS_gas_kgco2e2018_pdw
## t = -5.5938, df = 79, p-value = 3.081e-07
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.6727527 -0.3556615
## sample estimates:
##        cor 
## -0.5326464

Figure 5.4 uses the same plotting method to show emissions per dwelling due to electricity use.

## Per dwelling T CO2e - elec emissions
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Scatter of LSOA level elec per dwelling emissions against IMD score - who emits?

Figure 5.4: Scatter of LSOA level elec per dwelling emissions against IMD score - who emits?

## Correlation test (pwcorr)
## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$IMDScore and selectedLsoasDT$CREDS_elec_kgco2e2018_pdw
## t = 0.78479, df = 79, p-value = 0.4349
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1329488  0.3005317
## sample estimates:
##        cor 
## 0.08795387

Figure 5.5 uses the same plotting method to show emissions per dwelling due to other energy use. This should be higher for off-gas areas which tend to be rural areas so we also present analysis by the LSOA’s urban/rural classification.

## Per dwelling T CO2e - elec emissions
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Scatter of LSOA level other energy per dwelling emissions against IMD score - who emits?

Figure 5.5: Scatter of LSOA level other energy per dwelling emissions against IMD score - who emits?

## Correlation test (pwcorr)
## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$IMDScore and selectedLsoasDT$CREDS_elec_kgco2e2018_pdw
## t = 0.78479, df = 79, p-value = 0.4349
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.1329488  0.3005317
## sample estimates:
##        cor 
## 0.08795387
Table 5.5: Summary of per dwelling energy emissions for LSOAs by urban/rural code

RUC11

mean_gas_kgco2e

mean_elec_kgco2e

mean_other_energy_kgco2e

Rural town and fringe

2,435.3

1,292.6

373.0

Rural village and dispersed

2,228.3

1,774.2

976.7

Urban city and town

2,711.2

1,100.8

55.6

Urban major conurbation

2,421.5

1,205.4

25.9

Check whether all measured energy emissions combined (gas & electricity) correlate with all emissions (in this data).

## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$CREDS_total_kgco2e_pdw and selectedLsoasDT$CREDS_measuredHomeEnergy_kgco2e2018_pdw
## t = 8.3815, df = 79, p-value = 1.571e-12
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.5501273 0.7865839
## sample estimates:
##       cor 
## 0.6860664
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

Do we see strong correlations? If so in theory we could (currently) use measured energy emissions as a proxy for total emissions.

Repeat for all home energy - includes estimates of emissions from oil etc

## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$CREDS_total_kgco2e_pdw and selectedLsoasDT$CREDS_allHomeEnergy_kgco2e2018_pdw
## t = 8.342, df = 79, p-value = 1.877e-12
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  0.5478642 0.7853455
## sample estimates:
##       cor 
## 0.6843479
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

How does the correlation look now?

Figure 5.6: Annual mean T CO2e due to energy use per dwelling (LSOAs)

5.1.3 Transport

We don’t expect to use this data as it is already taxed in a way that relates to emissions (?)

Figure 5.7 uses the same plotting method to show emissions per dwelling due to van use. Again, we present analysis by the LSOA’s urban/rural classification.

## Per dwelling T CO2e - car emissions
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Scatter of LSOA level car use per dwelling emissions against IMD score

Figure 5.7: Scatter of LSOA level car use per dwelling emissions against IMD score

## Correlation with IMD score (pwcorr)
## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$IMDScore and selectedLsoasDT$CREDS_car_kgco2e2018_pdw
## t = -2.0208, df = 79, p-value = 0.04669
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.419731000 -0.003520997
## sample estimates:
##     cor 
## -0.2217
Table 5.6: Summary of per dwelling transport emissions for LSOAs by urban/rural code

RUC11

mean_car_kgco2e

mean_van_kgco2e

Rural town and fringe

3,370.6

515.5

Rural village and dispersed

3,877.9

685.5

Urban city and town

2,806.8

353.1

Urban major conurbation

3,004.9

379.9

Figure 5.8 uses the same plotting method to show emissions per dwelling due to van use.

## Per dwelling T CO2e - van emissions
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Scatter of LSOA level van use per dwelling emissions against IMD score

Figure 5.8: Scatter of LSOA level van use per dwelling emissions against IMD score

## Correlation with IMD score (pwcorr)
## 
##  Pearson's product-moment correlation
## 
## data:  selectedLsoasDT$IMDScore and selectedLsoasDT$CREDS_van_kgco2e2018_pdw
## t = 1.2808, df = 79, p-value = 0.204
## alternative hypothesis: true correlation is not equal to 0
## 95 percent confidence interval:
##  -0.0781545  0.3500759
## sample estimates:
##       cor 
## 0.1426289

Figure 5.9: Mean annual T CO2e per dwelling due to transport (LSOAs)

6 Estimating the annual emissions levy

Case studies:

  • Annual carbon tax
  • Half-hourly (real time) carbon tax (not implemented) - this would only affect electricity

BEIS/ETC Carbon ‘price’

EU carbon ‘price’

BEIS Carbon ‘Value’ https://www.gov.uk/government/publications/valuing-greenhouse-gas-emissions-in-policy-appraisal/valuation-of-greenhouse-gas-emissions-for-policy-appraisal-and-evaluation#annex-1-carbon-values-in-2020-prices-per-tonne-of-co2

  • based on a Marginal Abatement Cost (MAC)
  • 2021:
    • Low: £122/T
    • Central: £245/T <- use the central value for now
    • High: £367/T

Scenario 1: apply the central value Scenario 2: apply the low/central/high as a rising block tariff for each emissions source. Set threhsolds to 33% and 66% (in absence of any other guidance!)

6.1 Scenario 1: Central carbon cost

Table 6.1 below shows the overall £ GBP emissions levy total for the case study area in £M under Scenario 1.

Table 6.1: Total £m GBP levy by source (Scenario 1)

nLSOAs

All emissions

Gas

Electricity

81

344.4

37.9

18.5

Table 6.1: Total £m GBP levy by source, ordered by highest emitting ward (Scenario 1)

WD20NM

nLSOAs

All emissions

Gas

Electricity

Bishop's Stortford South

6

25.4

2.6

1.0

Hertford Bengeo

5

22.7

2.4

0.9

Hertford Castle

5

22.4

2.6

1.4

Bishop's Stortford Central

6

21.9

2.5

1.0

Bishop's Stortford All Saints

5

21.7

2.5

1.0

Sawbridgeworth

5

20.9

2.8

1.1

Bishop's Stortford Silverleys

4

19.0

1.9

0.6

Hertford Kingsmead

3

13.9

1.5

0.6

Hertford Sele

4

13.6

1.5

0.6

Bishop's Stortford Meads

3

13.4

2.0

0.9

Buntingford

3

12.7

1.9

0.8

Ware Christchurch

3

12.1

1.4

0.7

Ware St Mary's

3

11.8

1.3

0.6

Ware Trinity

3

10.8

1.4

0.5

Thundridge & Standon

2

8.8

0.7

0.5

Much Hadham

2

8.5

0.7

0.5

Ware Chadwell

2

8.3

0.8

0.5

Great Amwell

2

8.2

0.9

0.4

Hertford Heath

2

7.9

0.7

0.4

Hunsdon

2

7.7

0.8

0.5

Hertford Rural South

1

7.1

0.5

0.5

Stanstead Abbots

2

7.0

0.8

0.4

Little Hadham

1

6.0

0.3

0.5

Hertford Rural North

1

5.9

0.4

0.5

Datchworth & Aston

1

5.9

0.8

0.4

Puckeridge

1

4.6

0.7

0.4

Mundens and Cottered

1

4.5

0.1

0.3

Watton-at-Stone

1

4.5

0.6

0.3

Walkern

1

3.8

0.4

0.2

Braughing

1

3.6

0.3

0.2

Table 6.1: Total £m GBP levy by source, ordered by highest emitting ward (Scenario 1)

WD20NM

nLSOAs

All emissions

Gas

Electricity

Bishop's Stortford South

6

25.4

2.6

1.0

Hertford Bengeo

5

22.7

2.4

0.9

Hertford Castle

5

22.4

2.6

1.4

Bishop's Stortford Central

6

21.9

2.5

1.0

Bishop's Stortford All Saints

5

21.7

2.5

1.0

Sawbridgeworth

5

20.9

2.8

1.1

Bishop's Stortford Silverleys

4

19.0

1.9

0.6

Hertford Kingsmead

3

13.9

1.5

0.6

Hertford Sele

4

13.6

1.5

0.6

Bishop's Stortford Meads

3

13.4

2.0

0.9

Buntingford

3

12.7

1.9

0.8

Ware Christchurch

3

12.1

1.4

0.7

Ware St Mary's

3

11.8

1.3

0.6

Ware Trinity

3

10.8

1.4

0.5

Thundridge & Standon

2

8.8

0.7

0.5

Much Hadham

2

8.5

0.7

0.5

Ware Chadwell

2

8.3

0.8

0.5

Great Amwell

2

8.2

0.9

0.4

Hertford Heath

2

7.9

0.7

0.4

Hunsdon

2

7.7

0.8

0.5

Hertford Rural South

1

7.1

0.5

0.5

Stanstead Abbots

2

7.0

0.8

0.4

Little Hadham

1

6.0

0.3

0.5

Hertford Rural North

1

5.9

0.4

0.5

Datchworth & Aston

1

5.9

0.8

0.4

Puckeridge

1

4.6

0.7

0.4

Mundens and Cottered

1

4.5

0.1

0.3

Watton-at-Stone

1

4.5

0.6

0.3

Walkern

1

3.8

0.4

0.2

Braughing

1

3.6

0.3

0.2

Emissions level due to total emissions, gas & electricity use by LSOA (Scenario 1)

Figure 6.1: Emissions level due to total emissions, gas & electricity use by LSOA (Scenario 1)

The table below shows the mean per dwelling value rounded to the nearest £10.

Table 6.2: Mean £ GBP levy per dwelling by source (Scenario 1)

All_emissions

Gas

Electricity

Gas + Electricity

5,993.6

636.8

301.5

938.3

Figure ?? shows the total £k per LSOA and £ per dwelling revenue using BEIS central carbon price plotted against IMD score. The tables show the LSOAs with the highest and lowest values.

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA revenue using BEIS central carbon price

Figure 6.2: £k per LSOA revenue using BEIS central carbon price

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA revenue using BEIS central carbon price

Figure 6.3: £k per LSOA revenue using BEIS central carbon price

## Saving 7 x 5 in image
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##    2429    4325    6062    5994    6971   11819
Table 6.3: Highest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_allEmissions_levy

msoa_tot_annual_income_2018

E01023452

East Hertfordshire 003A

Bishop's Stortford Silverleys

48.2

11,818.8

56,600

E01023490

East Hertfordshire 015D

Hertford Sele

47.2

11,557.1

47,300

E01023468

East Hertfordshire 018A

Great Amwell

43.4

10,642.0

48,700

E01023461

East Hertfordshire 008F

Bishop's Stortford South

43.4

10,631.2

56,500

E01023471

East Hertfordshire 014B

Hertford Bengeo

41.9

10,269.0

56,600

E01023481

East Hertfordshire 018D

Hertford Heath

40.2

9,849.3

48,700

Table 6.3: Lowest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_allEmissions_levy

msoa_tot_annual_income_2018

E01023487

East Hertfordshire 015A

Hertford Sele

9.9

2,429.3

47,300

E01023450

East Hertfordshire 004B

Bishop's Stortford Meads

10.8

2,643.5

50,600

E01023499

East Hertfordshire 011A

Sawbridgeworth

10.8

2,654.3

49,800

E01023488

East Hertfordshire 015B

Hertford Sele

11.0

2,687.7

47,300

E01023478

East Hertfordshire 017D

Hertford Castle

11.5

2,810.7

53,000

E01023446

East Hertfordshire 007D

Bishop's Stortford Central

11.5

2,818.6

47,300

Figure ?? repeats the analysis but just for gas.

Anything unusual?

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA incurred via gas using BEIS central carbon price

Figure 6.4: £k per LSOA incurred via gas using BEIS central carbon price

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA incurred via gas using BEIS central carbon price

Figure 6.5: £k per LSOA incurred via gas using BEIS central carbon price

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   185.4   542.6   636.0   636.8   735.2  1068.0
Table 6.4: Highest emitting LSOAs by £ GBP gas levy

LSOA11CD

LSOA01NM

WD18NM

gas_T_CO2e_pdw

GBP_gas_levy_perdw

msoa_tot_annual_income_2018

E01023500

East Hertfordshire 011B

Sawbridgeworth

4.4

1,068.0

49,800

E01023474

East Hertfordshire 014E

Hertford Bengeo

3.9

952.3

56,600

E01023453

East Hertfordshire 003B

Bishop's Stortford Silverleys

3.9

948.5

56,600

E01023440

East Hertfordshire 005C

Bishop's Stortford All Saints

3.8

941.8

54,300

E01023460

East Hertfordshire 008E

Bishop's Stortford South

3.8

929.2

56,500

E01023454

East Hertfordshire 003C

Bishop's Stortford Silverleys

3.8

928.9

56,600

Table 6.4: Lowest emitting LSOAs by £ GBP gas levy

LSOA11CD

LSOA01NM

WD18NM

gas_T_CO2e_pdw

GBP_gas_levy_perdw

msoa_tot_annual_income_2018

E01023485

East Hertfordshire 010A

Hertford Rural North

1.7

406.1

56,700

E01023478

East Hertfordshire 017D

Hertford Castle

1.5

362.8

53,000

E01023493

East Hertfordshire 003F

Little Hadham

1.4

348.6

56,600

E01023496

East Hertfordshire 009C

Much Hadham

1.2

304.6

48,300

E01023511

East Hertfordshire 013B

Ware Chadwell

1.1

274.4

50,300

E01023494

East Hertfordshire 002B

Mundens and Cottered

0.8

185.4

53,300

Figure ?? repeats the analysis for electricity.

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA incurred via electricity using BEIS central carbon price

Figure 6.6: £k per LSOA incurred via electricity using BEIS central carbon price

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA incurred via electricity using BEIS central carbon price

Figure 6.7: £k per LSOA incurred via electricity using BEIS central carbon price

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   193.4   253.4   275.3   301.5   335.2   565.8
Table 6.5: Highest emitting LSOAs by £ GBP electricity levy

LSOA11CD

LSOA01NM

WD18NM

elec_T_CO2e_pdw

GBP_elec_levy_perdw

msoa_tot_annual_income_2018

E01023485

East Hertfordshire 010A

Hertford Rural North

2.3

565.8

56,700

E01023494

East Hertfordshire 002B

Mundens and Cottered

2.2

537.0

53,300

E01023493

East Hertfordshire 003F

Little Hadham

2.1

508.6

56,600

E01023486

East Hertfordshire 010B

Hertford Rural South

2.0

482.2

56,700

E01023507

East Hertfordshire 006C

Thundridge & Standon

1.9

453.9

48,700

E01023496

East Hertfordshire 009C

Much Hadham

1.8

440.1

48,300

Table 6.5: Lowest emitting LSOAs by £ GBP electricity levy

LSOA11CD

LSOA01NM

WD18NM

elec_T_CO2e_pdw

GBP_elec_levy_perdw

msoa_tot_annual_income_2018

E01023470

East Hertfordshire 014A

Hertford Bengeo

0.9

228.8

56,600

E01023447

East Hertfordshire 007E

Bishop's Stortford Central

0.9

220.5

47,300

E01023475

East Hertfordshire 017A

Hertford Castle

0.8

207.9

53,000

E01023489

East Hertfordshire 015C

Hertford Sele

0.8

205.8

47,300

E01023446

East Hertfordshire 007D

Bishop's Stortford Central

0.8

203.7

47,300

E01023487

East Hertfordshire 015A

Hertford Sele

0.8

193.4

47,300

Figure ?? shows the same analysis for measured energy (elec + gas)

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA incurred via electricity and gas using BEIS central carbon price

Figure 6.8: £k per LSOA incurred via electricity and gas using BEIS central carbon price

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
£k per LSOA incurred via electricity and gas using BEIS central carbon price

Figure 6.9: £k per LSOA incurred via electricity and gas using BEIS central carbon price

## Saving 7 x 5 in image
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   584.6   839.5   920.8   938.3  1023.0  1442.2

6.2 Scenario 2: Rising block tariff

Applied to per dwelling values (not LSOA total) - may be methodologically dubious?

Cut at 25%, 50% - so any emissions over 50% get high carbon cost

## Cuts for total per dw
##        0%       25%       50%       75%      100% 
##  9.915705 17.653546 24.741085 28.453846 48.240000
## `stat_bin()` using `bins = 30`. Pick better value with `binwidth`.

Table: (#tab:estimateAnnualLevyScenario2Total)Data summary

Name …[]
Number of rows 81
Number of columns 3
Key NULL
_______________________
Column type frequency:
numeric 3
________________________
Group variables None

Variable type: numeric

skim_variable n_missing complete_rate mean sd p0 p25 p50 p75 p100 hist
V1 0 1 24.46 9.04 9.92 17.65 24.74 28.45 48.24 ▆▇▇▂▂
beis_GBPtotal_sc2_perdw 0 1 4332.52 2731.78 1209.72 2153.73 3879.23 5252.76 12514.28 ▇▆▂▁▁
beis_GBPtotal_sc2 0 1 2944405.73 1376380.59 901238.40 1974391.43 2805946.21 3739013.30 6059048.95 ▆▇▆▂▂

Gas cuts:

## Cuts for gas per dw
##        0%       25%       50%       75%      100% 
## 0.7566615 2.2148481 2.5959184 3.0008683 4.3593750
## [1] 23.96836
## Cuts for elec per dw
##        0%       10%       20%       30%       40%       50%       60%       70%       80% 
## 0.7895839 0.9567590 1.0135132 1.0567619 1.1002625 1.1238281 1.1745586 1.2866934 1.4216188 
##       90%      100% 
## 1.6160106 2.3095535
## [1] 13.05225
Table 6.7: Highest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_allEmissions_levy

E01023452

East Hertfordshire 003A

Bishop's Stortford Silverleys

48.2

12,514.3

E01023490

East Hertfordshire 015D

Hertford Sele

47.2

12,122.3

E01023468

East Hertfordshire 018A

Great Amwell

43.4

10,751.6

E01023461

East Hertfordshire 008F

Bishop's Stortford South

43.4

10,735.4

E01023471

East Hertfordshire 014B

Hertford Bengeo

41.9

10,192.7

E01023481

East Hertfordshire 018D

Hertford Heath

40.2

9,564.0

Table 6.7: Lowest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_allEmissions_levy

E01023487

East Hertfordshire 015A

Hertford Sele

9.9

1,209.7

E01023450

East Hertfordshire 004B

Bishop's Stortford Meads

10.8

1,316.4

E01023499

East Hertfordshire 011A

Sawbridgeworth

10.8

1,321.7

E01023488

East Hertfordshire 015B

Hertford Sele

11.0

1,338.4

E01023478

East Hertfordshire 017D

Hertford Castle

11.5

1,399.6

E01023446

East Hertfordshire 007D

Bishop's Stortford Central

11.5

1,403.5

Table 6.7: Highest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_gasEmissions_levy

E01023500

East Hertfordshire 011B

Sawbridgeworth

4.4

1,010.8

E01023474

East Hertfordshire 014E

Hertford Bengeo

3.9

837.4

E01023453

East Hertfordshire 003B

Bishop's Stortford Silverleys

3.9

831.7

E01023440

East Hertfordshire 005C

Bishop's Stortford All Saints

3.8

821.7

E01023460

East Hertfordshire 008E

Bishop's Stortford South

3.8

802.8

E01023454

East Hertfordshire 003C

Bishop's Stortford Silverleys

3.8

802.3

Table 6.7: Lowest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_gasEmissions_levy

E01023494

East Hertfordshire 002B

Mundens and Cottered

0.8

92.3

E01023511

East Hertfordshire 013B

Ware Chadwell

1.1

136.6

E01023496

East Hertfordshire 009C

Much Hadham

1.2

151.7

E01023493

East Hertfordshire 003F

Little Hadham

1.4

173.6

E01023478

East Hertfordshire 017D

Hertford Castle

1.5

180.6

E01023485

East Hertfordshire 010A

Hertford Rural North

1.7

202.2

Table 6.7: Highest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_elecEmissions_levy

E01023485

East Hertfordshire 010A

Hertford Rural North

2.3

606.3

E01023494

East Hertfordshire 002B

Mundens and Cottered

2.2

563.1

E01023493

East Hertfordshire 003F

Little Hadham

2.1

520.6

E01023486

East Hertfordshire 010B

Hertford Rural South

2.0

481.0

E01023507

East Hertfordshire 006C

Thundridge & Standon

1.9

438.6

E01023496

East Hertfordshire 009C

Much Hadham

1.8

417.9

Table 6.7: Lowest levy LSOAs

LSOA11CD

LSOA01NM

WD18NM

T_CO2e_pdw

GBP_elecEmissions_levy

E01023487

East Hertfordshire 015A

Hertford Sele

0.8

96.3

E01023446

East Hertfordshire 007D

Bishop's Stortford Central

0.8

101.4

E01023489

East Hertfordshire 015C

Hertford Sele

0.8

102.5

E01023475

East Hertfordshire 017A

Hertford Castle

0.8

103.5

E01023447

East Hertfordshire 007E

Bishop's Stortford Central

0.9

109.8

E01023470

East Hertfordshire 014A

Hertford Bengeo

0.9

113.9

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

Table 6.8: Total £m GBP levy by source (Scenario 2)

nLSOAs

All emissions

Gas

Electricity

81

238.5

24.0

13.1

Table 6.8: Total £m GBP levy by source, ordered by highest emitting ward (Scenario 2)

WD20NM

nLSOAs

All emissions

Gas

Electricity

Bishop's Stortford South

6

21.0

1.9

0.6

Hertford Bengeo

5

17.9

1.7

0.6

Bishop's Stortford Silverleys

4

16.8

1.5

0.4

Bishop's Stortford Central

6

14.8

1.4

0.5

Sawbridgeworth

5

14.7

2.0

0.7

Bishop's Stortford All Saints

5

14.3

1.8

0.6

Hertford Castle

5

14.1

1.4

0.9

Hertford Sele

4

10.6

0.9

0.3

Hertford Kingsmead

3

9.3

0.9

0.3

Bishop's Stortford Meads

3

8.1

1.4

0.6

Great Amwell

2

7.1

0.6

0.3

Buntingford

3

7.0

1.1

0.5

Ware St Mary's

3

6.9

0.6

0.4

Ware Christchurch

3

6.7

0.8

0.4

Thundridge & Standon

2

6.6

0.4

0.5

Ware Trinity

3

6.4

0.8

0.3

Much Hadham

2

6.4

0.5

0.5

Hertford Heath

2

6.2

0.5

0.3

Ware Chadwell

2

5.3

0.5

0.4

Hertford Rural South

1

5.0

0.3

0.5

Hunsdon

2

4.6

0.4

0.5

Stanstead Abbots

2

4.3

0.4

0.2

Hertford Rural North

1

3.9

0.2

0.6

Little Hadham

1

3.9

0.2

0.5

Datchworth & Aston

1

3.7

0.7

0.3

Mundens and Cottered

1

3.4

0.1

0.4

Walkern

1

2.6

0.2

0.2

Braughing

1

2.4

0.2

0.2

Puckeridge

1

2.3

0.3

0.3

Watton-at-Stone

1

2.2

0.3

0.2

Total levy due to all emissions, gas & electricity use by LSOA (Scenario 2)

Figure 6.10: Total levy due to all emissions, gas & electricity use by LSOA (Scenario 2)

The table below shows the mean per dwelling value rounded to the nearest £10.

Table 6.9: Mean £ GBP levy per dwelling by source (Scenario 2)

All_emissions

Gas

Electricity

1,078.4

411.2

213.0

6.3 Compare scenarios

Figure 6.11 compares the % £ levy under each scenario for all consumption contributed by LSOAs in each IMD decile.

Table 6.10: Compare totals for Scenario 1 & Scenario 2 (£m)

nLSOAs

sum_total_sc1

sum_total_sc2

81

344.4

238.5

Comparing £ levy under each scenario by IMD decile - all consumption emissions

Figure 6.11: Comparing £ levy under each scenario by IMD decile - all consumption emissions

## Saving 7 x 5 in image

Figure 6.12 compares the £ levy under each scenario for all consumption.

## Scale for colour is already present.
## Adding another scale for colour, which will replace the existing scale.
Comparing £ levy under each scenario - all consumption emissions

Figure 6.12: Comparing £ levy under each scenario - all consumption emissions

7 Estimate retofit costs

  • from A-E <- £13,300
  • from F-G <- £26,800

Source: English Housing Survey 2018 Energy Report

Model excludes EPC A, B & C (assumes no need to upgrade)

Adding these back in would increase the cost… obvs

7.1 Impute EPC counts

There are likely to be a range of biases in the EPC sample such as:

  • required for new builds
  • required for rental
  • required for sale
  • required after some retrofit

Check distributions of EPCs against known building stock:

  • Census dwelling type vs EPC dwelling type
  • Census dwelling age (?) vs EPC building age
  • Census main heating type ? vs EPC “mainfuel_mainsgas” “mainfuel_electric” “mainfuel_oil” “mainfuel_coal” “mainfuel_lpg” “mainfuel_biomass”

##  pc_missing_epcs
##  Min.   :14.36  
##  1st Qu.:42.16  
##  Median :48.04  
##  Mean   :45.95  
##  3rd Qu.:53.00  
##  Max.   :57.34

These plots suggest:

  • Lots of detached houses missing - especially in lower deprivation areas
  • Lots of semis missing
  • Terrace houses better represented
  • Flats fairly well represented
  • Overall compared to the number of electricity meters, we have about 40% of the EPCs we should have

In order to estimate the LSOA level retrofit costs, we need to impute the EPC counts in each LSOA. We do this using the number of electricity meters as the presumed number of dwellings and the observed % of EPCs in each band for all dwellings with EPCs which is provided by the CREDS data. We simply upweight the EPC counts proportionally so that the total matches the number of electricity meters.

Note that this assumes that if we had EPCs for all dwellings then the % in each band in each LSOA would stay the same. This is quite a bold assumption given what we have seen above but there is currently no way to re-weight them (e.g. to fit the dwelling type counts) with the data we have - we do not know the of EPCs of each Band that come from different dwelling types.

Note to self: this could be done by re-weighting the case level EPC data using Census outcomes at LSOA level - e.g. via IPF

Note that the EPC database is continuously updated so more recent upgrades will not be captured in the data used for this analysis. This means the total retrofit costs are likely to be an over-estimate. The extent of this over-estimate would require the use of an updated (current) EPC data extract and is left for future work.

## N EPCs
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   214.0   291.0   358.0   421.5   504.0  1360.0
## N elec meters
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##     469     606     690     758     885    1676
## [1] 100
## [1] 100

variable

nDwellings

pc

i.nDwellings

i.pc

nDwellings_A

44.84

0.07

12,737.80

0.17

nDwellings_B

9,160.00

15.18

778,932.47

10.59

nDwellings_C

17,611.83

29.18

2,208,676.15

30.04

nDwellings_D

22,906.20

37.95

2,939,176.13

39.98

nDwellings_E

8,938.41

14.81

881,518.85

11.99

nDwellings_F

1,172.00

1.94

502,229.00

6.83

nDwellings_G

519.01

0.86

28,665.51

0.39

Correlation between high % EPC F/G or A/B and deprivation?

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

Now we need to convert the % to dwellings using the number of electricity meters (see above).

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

7.2 Estimate costs - sheep dip D to G

Upgrade all D-G:

  • upgrading A-E to ‘acceptable standard’: £13,300
  • upgrading F-G to ‘acceptable standard’: £26,800

Source: estimates of the cost of retrofitting dwellings in different EPC bands ‘to an appropriate standard’ published in a recent English Housing Survey technical report (MHCLG 2018) (see https://journal-buildingscities.org/articles/10.5334/bc.279#T2)

The model does not upgrade A-C as it was assumed that dwellings in the least energy-efficient EPC bands (D–G) would be prioritised and that EPC bands A–C would not be retrofitted since they would already conform to the current UK policy objective of:

“as many homes as possible to achieve EPC band C by 2035 where cost-effective, practical and affordable, and to ensure as many fuel poor homes as reasonably practicable achieve a band C rating by the end of 2030.” https://www.gov.uk/government/news/plan-to-drive-down-the-cost-of-clean-heat

Table 7.1 reports total retofit costs.

## To retrofit D-E (£m)
## [1] 423.5333
## Number of dwellings: 31845
## To retrofit F-G (£m)
## [1] 73.10761
## Number of dwellings: 2728
## To retrofit D-G (£m)
## [1] 496.6409
## To retrofit D-G (mean per dwelling)
## [1] 14280.66
Table 7.1: Retrofit cost totals (£m GBP):

Mean total £m per LSOA

Total £m

6.1

496.6

Table 7.1: Retrofit cost totals by ward (£m GBP):

Ward

Mean total £m per LSOA

Total £m

Sawbridgeworth

6.7

33.7

Hertford Castle

6.2

30.8

Hertford Bengeo

6.0

30.1

Bishop's Stortford South

4.9

29.5

Bishop's Stortford Central

4.9

29.4

Bishop's Stortford All Saints

5.6

28.1

Bishop's Stortford Meads

7.5

22.4

Ware Christchurch

7.0

20.9

Ware St Mary's

6.9

20.8

Ware Trinity

6.4

19.1

Buntingford

6.0

17.9

Hertford Kingsmead

5.9

17.6

Hertford Sele

4.4

17.5

Bishop's Stortford Silverleys

4.2

16.7

Hertford Rural North

12.9

12.9

Little Hadham

12.8

12.8

Thundridge & Standon

6.4

12.8

Much Hadham

6.3

12.5

Hertford Rural South

12.0

12.0

Hunsdon

5.9

11.8

Stanstead Abbots

5.7

11.4

Hertford Heath

5.4

10.8

Puckeridge

10.5

10.5

Datchworth & Aston

10.4

10.4

Ware Chadwell

5.0

9.9

Great Amwell

4.4

8.7

Watton-at-Stone

8.4

8.4

Mundens and Cottered

6.1

6.1

Braughing

5.8

5.8

Walkern

5.3

5.3

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

Figure 7.1 shows the LSOA level retofit costs per dwelling by IMD decile.

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
LSOA level retofit costs per dwelling by IMD score

Figure 7.1: LSOA level retofit costs per dwelling by IMD score

## Saving 7 x 5 in image

Figure 7.2: Total estimated retorfit costs per LSOA

7.2.1 Estimate costs - sheep dip D to G but upgrades in 20% poorest areas funded by HM Gov (somehow)

Set aside areas in IMD decile 1 & 2? Assume costs met by HM Treasury or elsewhere e.g. energy levy etc

to do

8 Compare levy with costs

8.1 Scenario 1

Totals

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

Repeat per dwelling

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

8.2 Scenario 2

Totals

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

Repeat per dwelling

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'

9 Years to pay…

9.1 Scenario 1

Figure 9.1 shows years to pay under Scenario 1 (all emissions)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   1.132   2.101   2.461   2.742   3.350   5.546
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Years to pay under Scenario 1 (all em issions)

Figure 9.1: Years to pay under Scenario 1 (all em issions)

## Saving 7 x 5 in image
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
## Median years: 2.46

Figure 9.2 maps the distribution of years to pay for all LSOAs in the case study.

Figure 9.2: Years to pay per LSOA (Scenbario 1, all emissions levy, no redistribution)

Figure 9.3 shows years to pay under Scenario 1 (energy emissions)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   9.621  13.682  15.714  15.746  17.136  24.180
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Years to pay under Scenario 1 (energy emissions)

Figure 9.3: Years to pay under Scenario 1 (energy emissions)

## Median years: 15.71

Figure 9.4 shows those LSOAs which would be in surplus (above the y = 0 line) after Year 1.

Surplus £ after year 1 for each LSOA (all emissions levy)

Figure 9.4: Surplus £ after year 1 for each LSOA (all emissions levy)

Figure 9.5 shows the same distribution but for the energy emissions levy. Clearly there is a lower level of surplus due to the reduced levy value at LSOA level.

Surplus £ after year 1 for each LSOA (energy levy)

Figure 9.5: Surplus £ after year 1 for each LSOA (energy levy)

What happens in Year 2 totally depends on the rate of upgrades… given the supply chain & capacity issues it’s likely that the levy would build up a substantial ‘headroom’ that could then be spent over time…

Figure 9.6 shows the year 1 outcome at ward level if levy is shared equally across all LSOAs (all emissions levy - orange dots).

Year 1 outcome if levy is shared equally (all emissions levy)

Figure 9.6: Year 1 outcome if levy is shared equally (all emissions levy)

## Saving 7 x 5 in image
Table 9.1: Longest years to pay under Scenario 1 (sorted by total retrofit cost)

LSOA

LSOA name

Ward

Retrofit cost £m

Years to pay (all emissions levy)

Years to pay (energy emissions levy)

% EPC D

% EPC E

% EPC F

% EPC G

E01023485

East Hertfordshire 010A

Hertford Rural North

12.9

2.5

16.0

41.2

30.2

12.4

2.3

E01023493

East Hertfordshire 003F

Little Hadham

12.8

2.8

19.5

31.9

25.9

17.1

2.1

E01023486

East Hertfordshire 010B

Hertford Rural South

12.0

2.5

16.9

33.3

21.6

9.7

3.6

E01023498

East Hertfordshire 006A

Puckeridge

10.5

3.9

16.5

38.4

15.9

3.4

0.9

E01023467

East Hertfordshire 002A

Datchworth & Aston

10.4

2.5

12.1

45.9

18.8

5.8

1.5

E01023516

East Hertfordshire 013D

Ware St Mary's

8.5

2.7

16.4

53.4

13.8

2.3

0.3

E01033609

East Hertfordshire 010E

Watton-at-Stone

8.4

3.7

16.9

36.4

9.9

4.1

1.1

E01023451

East Hertfordshire 004C

Bishop's Stortford Meads

8.4

2.9

13.6

37.2

17.6

3.1

0.2

E01023446

East Hertfordshire 007D

Bishop's Stortford Central

8.2

5.1

23.2

24.0

15.2

2.6

0.9

E01023514

East Hertfordshire 012C

Ware Christchurch

8.0

4.4

21.1

26.2

16.6

4.1

1.6

9.2 Scenario 2

Figure 9.7 shows years to pay under Scenario 2 (all emissions)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   1.069   2.786   3.799   4.689   6.727  11.137
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Years to pay under Scenario 2 (all emissions)

Figure 9.7: Years to pay under Scenario 2 (all emissions)

## Median years: 3.8

Figure 9.8 shows years to pay under Scenario 2 (energy emissions)

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   10.43   19.59   24.91   25.45   30.12   46.50
## `geom_smooth()` using method = 'loess' and formula = 'y ~ x'
Years to pay under Scenario 2 (energy emissions)

Figure 9.8: Years to pay under Scenario 2 (energy emissions)

Figure 9.9 shows those LSOAs which would be in surplus (above the y = 0 line) after Year 1.

Surplus £ after year 1 for each LSOA (Scenario 2, all emissions levy)

Figure 9.9: Surplus £ after year 1 for each LSOA (Scenario 2, all emissions levy)

Figure 9.10 shows the same distribution for the Scenario 2 energy levy.

Surplus £ after year 1 for each LSOA (Scenario 2, all emissions levy)

Figure 9.10: Surplus £ after year 1 for each LSOA (Scenario 2, all emissions levy)

What happens in Year 2 totally depends on the rate of upgrades…

Figure 9.11 shows the year 1 outcome at ward level if levy is shared equally across all LSOAs (all emissions levy - orange dots).

Year 1 outcome if levy is shared equally (Scenario 2, all emissions levy)

Figure 9.11: Year 1 outcome if levy is shared equally (Scenario 2, all emissions levy)

## Saving 7 x 5 in image
Table 9.2: Longest years to pay under Scenario 2 (sorted by total retrofit cost)

LSOA

LSOA name

Ward

Retrofit cost £m

Years to pay (all emissions levy)

Years to pay (energy emissions levy)

% EPC D

% EPC E

% EPC F

% EPC G

E01023485

East Hertfordshire 010A

Hertford Rural North

12.9

3.9

19.3

41.2

30.2

12.4

2.3

E01023493

East Hertfordshire 003F

Little Hadham

12.8

4.3

24.0

31.9

25.9

17.1

2.1

E01023486

East Hertfordshire 010B

Hertford Rural South

12.0

3.5

22.4

33.3

21.6

9.7

3.6

E01023498

East Hertfordshire 006A

Puckeridge

10.5

7.9

27.9

38.4

15.9

3.4

0.9

E01023467

East Hertfordshire 002A

Datchworth & Aston

10.4

3.9

14.8

45.9

18.8

5.8

1.5

E01023516

East Hertfordshire 013D

Ware St Mary's

8.5

4.7

29.7

53.4

13.8

2.3

0.3

E01033609

East Hertfordshire 010E

Watton-at-Stone

8.4

7.4

29.0

36.4

9.9

4.1

1.1

E01023451

East Hertfordshire 004C

Bishop's Stortford Meads

8.4

5.2

19.6

37.2

17.6

3.1

0.2

E01023446

East Hertfordshire 007D

Bishop's Stortford Central

8.2

10.3

46.5

24.0

15.2

2.6

0.9

E01023514

East Hertfordshire 012C

Ware Christchurch

8.0

8.7

38.4

26.2

16.6

4.1

1.6

9.3 Compare scenarios

Figure 9.12 compares pay-back times for the two scenarios - who does the rising block tariff help?

Comparing pay-back times across scenarios

Figure 9.12: Comparing pay-back times across scenarios

## Saving 7 x 5 in image

10 Summary: East Hertfordshire

Table 10.1 shows the summary overall levy value results for East Hertfordshire by Scenario.

Table 10.1: Levy total (Year 1 £m GBP): East Hertfordshire

Scenario

nLSOAs

All emissions levy

Gas emissions levy

Electricity emissions levy

Scenario 1

81

344.4

37.9

18.5

Scenario 2

81

238.5

24.0

13.1

Table 10.2 shows the summary overall retrofit cost results for East Hertfordshire.

Table 10.2: Retrofit cost totals (£m GBP): East Hertfordshire

Total retrofit cost £m

Mean total retrofit cost £m per LSOA

496.6

6.1

As a reminder:

  • Number of households (Census 2021): 53,505
  • Estimated number of domestic electricity meters (2018): 61,397
  • Estimated number of domestic gas meters (2018): 50,756

The following table (Table 10.3) summarises all results by ward.

Table 10.3: All results by ward (£m GBP, sorted by total retrofit cost):

WD20NM

nLSOAs

All emissions levy (£m, Scenario 1)

Gas emissions levy (£m, Scenario 1)

Electricity emissions levy (£m, Scenario 1)

All emissions levy (£m, Scenario 2)

Gas emissions levy (£m, (Scenario 2)

Electricity emissions levy (£m, Scenario 2)

Total retrofit cost (£m)

Mean total retrofit cost per LSOA (£m)

Sawbridgeworth

5

20.9

2.8

1.1

14.7

2.0

0.7

33.7

6.7

Hertford Castle

5

22.4

2.6

1.4

14.1

1.4

0.9

30.8

6.2

Hertford Bengeo

5

22.7

2.4

0.9

17.9

1.7

0.6

30.1

6.0

Bishop's Stortford South

6

25.4

2.6

1.0

21.0

1.9

0.6

29.5

4.9

Bishop's Stortford Central

6

21.9

2.5

1.0

14.8

1.4

0.5

29.4

4.9

Bishop's Stortford All Saints

5

21.7

2.5

1.0

14.3

1.8

0.6

28.1

5.6

Bishop's Stortford Meads

3

13.4

2.0

0.9

8.1

1.4

0.6

22.4

7.5

Ware Christchurch

3

12.1

1.4

0.7

6.7

0.8

0.4

20.9

7.0

Ware St Mary's

3

11.8

1.3

0.6

6.9

0.6

0.4

20.8

6.9

Ware Trinity

3

10.8

1.4

0.5

6.4

0.8

0.3

19.1

6.4

Buntingford

3

12.7

1.9

0.8

7.0

1.1

0.5

17.9

6.0

Hertford Kingsmead

3

13.9

1.5

0.6

9.3

0.9

0.3

17.6

5.9

Hertford Sele

4

13.6

1.5

0.6

10.6

0.9

0.3

17.5

4.4

Bishop's Stortford Silverleys

4

19.0

1.9

0.6

16.8

1.5

0.4

16.7

4.2

Hertford Rural North

1

5.9

0.4

0.5

3.9

0.2

0.6

12.9

12.9

Little Hadham

1

6.0

0.3

0.5

3.9

0.2

0.5

12.8

12.8

Thundridge & Standon

2

8.8

0.7

0.5

6.6

0.4

0.5

12.8

6.4

Much Hadham

2

8.5

0.7

0.5

6.4

0.5

0.5

12.5

6.3

Hertford Rural South

1

7.1

0.5

0.5

5.0

0.3

0.5

12.0

12.0

Hunsdon

2

7.7

0.8

0.5

4.6

0.4

0.5

11.8

5.9

Stanstead Abbots

2

7.0

0.8

0.4

4.3

0.4

0.2

11.4

5.7

Hertford Heath

2

7.9

0.7

0.4

6.2

0.5

0.3

10.8

5.4

Puckeridge

1

4.6

0.7

0.4

2.3

0.3

0.3

10.5

10.5

Datchworth & Aston

1

5.9

0.8

0.4

3.7

0.7

0.3

10.4

10.4

Ware Chadwell

2

8.3

0.8

0.5

5.3

0.5

0.4

9.9

5.0

Great Amwell

2

8.2

0.9

0.4

7.1

0.6

0.3

8.7

4.4

Watton-at-Stone

1

4.5

0.6

0.3

2.2

0.3

0.2

8.4

8.4

Mundens and Cottered

1

4.5

0.1

0.3

3.4

0.1

0.4

6.1

6.1

Braughing

1

3.6

0.3

0.2

2.4

0.2

0.2

5.8

5.8

Walkern

1

3.8

0.4

0.2

2.6

0.2

0.2

5.3

5.3

As ever, #YMMV.

11 R environment

11.1 R packages used

  • base R (R Core Team 2016)
  • bookdown (Xie 2016a)
  • data.table (Dowle et al. 2015)
  • flextable (Gohel 2021)
  • ggplot2 (Wickham 2009)
  • here (Müller 2017)
  • knitr (Xie 2016b)
  • leaflet (Cheng, Karambelkar, and Xie 2023)
  • lubridate (Grolemund and Wickham 2011)
  • rmarkdown (Allaire et al. 2018)
  • sf (Pebesma 2018)
  • skimr (Arino de la Rubia et al. 2017)
  • viridisLite (Garnier et al. 2023)

References

Allaire, JJ, Yihui Xie, Jonathan McPherson, Javier Luraschi, Kevin Ushey, Aron Atkins, Hadley Wickham, Joe Cheng, and Winston Chang. 2018. Rmarkdown: Dynamic Documents for r. https://CRAN.R-project.org/package=rmarkdown.
Arino de la Rubia, Eduardo, Hao Zhu, Shannon Ellis, Elin Waring, and Michael Quinn. 2017. Skimr: Skimr. https://github.com/ropenscilabs/skimr.
Cheng, Joe, Bhaskar Karambelkar, and Yihui Xie. 2023. Leaflet: Create Interactive Web Maps with the JavaScript ’Leaflet’ Library. https://CRAN.R-project.org/package=leaflet.
Dowle, M, A Srinivasan, T Short, S Lianoglou with contributions from R Saporta, and E Antonyan. 2015. Data.table: Extension of Data.frame. https://CRAN.R-project.org/package=data.table.
Garnier, Simon, Ross, Noam, Rudis, Robert, Camargo, et al. 2023. viridis(Lite) - Colorblind-Friendly Color Maps for r. https://doi.org/10.5281/zenodo.4678327.
Gohel, David. 2021. Flextable: Functions for Tabular Reporting. https://CRAN.R-project.org/package=flextable.
Grolemund, Garrett, and Hadley Wickham. 2011. “Dates and Times Made Easy with lubridate.” Journal of Statistical Software 40 (3): 1–25. http://www.jstatsoft.org/v40/i03/.
Müller, Kirill. 2017. Here: A Simpler Way to Find Your Files. https://CRAN.R-project.org/package=here.
Pebesma, Edzer. 2018. Simple Features for R: Standardized Support for Spatial Vector Data.” The R Journal. https://journal.r-project.org/archive/2018/RJ-2018-009/index.html.
R Core Team. 2016. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing. https://www.R-project.org/.
Wickham, Hadley. 2009. Ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York. http://ggplot2.org.
Xie, Yihui. 2016a. Bookdown: Authoring Books and Technical Documents with R Markdown. Boca Raton, Florida: Chapman; Hall/CRC. https://github.com/rstudio/bookdown.
———. 2016b. Knitr: A General-Purpose Package for Dynamic Report Generation in r. https://CRAN.R-project.org/package=knitr.